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Abstract 

Spinal cord injury (SCI) is a devastating event with a limited hope for recovery and represents an enormous public health 
issue. It is crucial to understand the disturbances in the metabolic network after SCI to identify injury mechanisms and 
opportunities for treatment intervention. Through plasma 1H-nuclear magnetic resonance (NMR) screening, we identified 
15 metabolites that made up an "Eigen-metabolome" capable of distinguishing rats with severe SCI from healthy control 
rats. Forty enzymes regulated these 15 metabolites in the metabolic network. We also found that 16 metabolites regulated 
by 130 enzymes in the metabolic network impacted neurobehavioral recovery. Using the Eigen-metabolome, we 
established a linear discrimination model to cluster rats with severe and mild SCI and control rats into separate groups and 
identify the interactive relationships between metabolic biomarkers in the global metabolic network. We identified 10 
clusters in the global metabolic network and defined them as distinct metabolic disturbance domains of SCI. Metabolic 
paths such as retinal, glycerophospholipid, arachidonic acid metabolism; NAD-NADPH conversion process, tyrosine 
metabolism, and cadaverine and putrescine metabolism were included. In summary, we presented a novel interdisciplinary 
method that integrates metabolomics and global metabolic network analysis to visualize metabolic network disturbances 
after SCI. Our study demonstrated the systems biological study paradigm that integration of 1H-NMR, metabolomics, and 
global metabolic network analysis is useful to visualize complex metabolic disturbances after severe SCI. Furthermore, our 
findings may provide a new quantitative injury severity evaluation model for clinical use. 
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Introduction 

Spinal cord injury (SCI) is a major public health challenge that 
often affects young adults. In the United States, 12,000 new SCI 
cases are annually reported, and $24 billion is annually spent on 
paralyzed SCI patients [1]. The incidence of SCI in China is 
surging because of the large population and rapid economic 
growth; in Beijing the incidence has reached 60/ 10 6 per year [2]. 
SCI is a complex injury that involves multiple pathological 
processes, particularly when the injury is severe, and triggers shock 
and organ dysfunction or death. Understanding the cellular and 
metabolic network malfunction during SCI is crucial for clinical 
monitoring and intervention. 

Human metabolism is a complex network with hundreds of 
cross-linked paths. During critical SCI illness, the metabolic 
network is dynamically disturbed at multiple points [3]. Classical 



research typically isolates a small part of this network to investigate 
the impact of pathophysiological molecular mechanisms on 
clinical outcome. In particular, researchers have examined 
metabolic disturbances such as skeletal muscle breakdown, insulin 
resistance, dyslipidemia, testosterone and growth hormone/IGF-I 
dysfunctions, low thyroxine syndrome, and deficiency of vitamin D 
and calcium with secondary hyperparathyroidism [3-14]. These 
complex metabolic disturbances appear and interact at different 
stages during the pathological process after SCI [3,4]. Therefore, 
an integrated approach that combines the biochemical/ molecular 
changes with network disturbances is the key to understanding 
SCI at the systems biology level and establishing an accurate 
quantitative model for monitoring SCI. 

An interdisciplinary method that includes high-throughput 
quantitative techniques and effective mathematical and visualiza- 
tion tools is necessary because metabolic network disturbances 
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after SCI are complex; hundreds of molecules and metabolic paths 
interact [15,16]. Furthermore, interdisciplinary methods present 
the opportunity to develop innovative clinical diagnosis and 
monitoring methods for severe injuries. Here we present a novel 
high-throughput method that integrated lH-nuclear magnetic 
resonance (NMR) metabolomic fingerprinting with global meta- 
bolic network analysis to formulate a quantitative SCI evaluation 
model. 

Materials and Methods 

1. Modeling Spinal Cord Injury in SD Rats 

Animal modeling procedures were performed according to an 
established protocol, as described previously [17]. In brief, adult 
male Sprague-Dawley (SD) rats that had a normal neurobehav- 
ioral function were included in the study. The Basso, Beattie, and 
Bresnahan (BBB) score scale was used to evaluate neurobehavioral 
function (BBB scores >20 were considered normal, BBB scores of 
enrolled animals see Supporting Information: table SI and table 
S3). Rats were randomly divided into three groups: control [n — 7), 
mild SCI (n — 10), and severe SCI (n = 25). All the rats were fed the 
standard formula feed from 7 days before the operation until 
sacrifice. All rat forage was formulated according to the 
formulation established by the State Bureau of Quality Technical 
Supervision of China (GB Standards No: 14924.3-2010). Animal 
use was strictly in accordance with the National Institutes of 
Health Guide for the Care and Use of Laboratory Animals. The 
Experimental Animal Committee of Sichuan Academy of Medical 
Sciences approved the research protocol (Approval Number: 10- 
1069). 

Rats assigned to the SCI group were anesthetized by intra- 
peritoneal injection of sodium pentobarbital (40 mg/kg). We tried 
our best to minimize animal suffering. Rats were observed for at 
least 3 days after the start of the experiment. For severe SCI, we 
established a 7-day observation subgroup. In the mild SCI group, 
three rats died during the surgical modeling procedure and one 
died 2 days after the surgery. In the severe SCI group, five rats 
died during the surgical modeling procedure and four died 1 and 2 
days after the surgery. Therefore, there were six animals in the 
mild SCI group, 16 in the severe injury group (seven rats were fed 
for 7 days after modeling), and seven in the blank control group. 
Plasma samples were collected for 1H-NMR metabolomic 
analysis. 

2. 1H-NMR Metabolomic Testing 

Plasma samples were centrifuged at 10,000 rpm for 10 min to 
remove particulate matter. Following this, 100 jiL of deuterium 
oxide (D20) was added to 400 jiL of the plasma samples in 5 mm 
Wilmad NMR tubes to lock the field frequency for the NMR 
spectrometer. 

A high-solution Bruker AVIII-DRX600MHz NMR spectrom- 
eter (600.13 MHz 1H-NMR observation frequency; Bruker 
Biospin Rheinstetten, Germany) was used to obtain 1H-NMR 
spectra from the plasma samples. The probe temperature was 300 
K. One 1D-NMR experiment was included. A standard protocol 
that included a one-dimensional pulse sequence and a Carr- 
Purcell-Meiboom-Gill (CPMG) pulse sequence was followed [34], 
All samples were subjected to a 64 ms spin-spin relaxation delay 
and a 2-s water suppression irradiation during the relaxation 
delay. The spectral width was 20 ppm, and 64 transients were 
collected into 64,000 data points. The CPMG pulse sequence was 
used to filter broad resonances from protein and lipids and to 
visualize the latent biomarkers of smaller molecules. 



3. Data processing and multivariate analysis 

Principle component analysis (PCA) and partial least-square 
(PLS) regression were performed on 1 H-NMR data to reduce data 
dimensionality and for multivariate analysis. We used Matlab 
R2012b for computing and mathematical modeling (The Math- 
works, Inc.; Natwick, MA, US). The variable importance for the 
projection (VIP) value was used to screen metabolites. The Human 
Metabolome Database (HMDB) was used to identify key 
metabolites related to enzymes and upstream genes. Metabolites 
and their enzymes were mapped, and the disturbed metabolic 
paths were located on a global metabolic map of the Kyoto 
Encyclopedia of Genes and Genomes (KEGG). 

A linear discriminant formula was established as a quantitative 
injury evaluation model. In this model, weighted PLS-DA values 
were used. The root mean square error (RMSE) was calculated to 
verify the validity of the equations according to injury severity. 

Metabolic disturbances after severe SCI were established 
according to two steps: 1) key metabolites were mapped onto the 
KEGG metabolic network (KEGG metabolic pathways PATH- 
WAY: mapO 1 1 00) and 2) the shortest path for each key metabolite 
(node) was calculated and outlined onto the metabolic network 
using the KEGG pathway API for the Matlab bioinformatics 
toolbox. 

Results 

1. Effects of SCI on the plasma metabolome 

Compared with the control group, mild and severe SCI 1H- 
NMR spectra were perturbed at 0-8 ppm (Figure 1). This 
phenomenon suggests that the profiles of low-molecular-weight 
plasma metabolites and lipoproteins were altered by SCI. Visual 
inspection of the spectra demonstrated that there were significant 
alterations in plasma lactate, lysine, glucose, choline, and lipids 
(Table 1 and Figure 1). 

a) PCA. PCA scores indicated that injured rat plasma 
metabolomes differed from controls. The first two principal 
components (PCs) explained 83% of the total variance and 
suggested that most metabolome differences between normal 
and SCI cases were because of injury (Figure 2 A). 

2. The PLS model: severe SCI Eigen-metabolomes 

PLS was used to investigate correlations between severity of SCI 
and the metabolome. The VIP indicator can describe correlations 
between the variable (X) and response (Y). We used VIP to identify 
metabolites correlated with severe SCI and named these 
metabolites as "Eigen-metabolome" of severe injury. Segmental 
integrations in 83, 90, 92, 151, 152, and 153 were found to have 
the highest VIP values and represented metabolites with 4.72, 
4.44, 4.4, 2, 1.98, and 1.96 1 H-NMR spectra chemical shift 
values, respectively. 

We identified 15 metabolites that made up an "Eigen- 
metabolome" capable of distinguishing rats with severe SCI from 
healthy control rats. We also found that 16 metabolites impacted 
neurobehavioral recovery (Table 1 and 2). 

3. Quantitative severe SCI discriminate model 

To develop a quantitative severe SCI evaluation model based 
on the Eigen-metabolome, PLS regression analysis was performed 
to establish a linear prediction model (equation 1). The accuracy 
ratio of PLS regression was 86.2%, R 2 = 0.284 (Figure 3F-G). X t 
represented the z-th ppm value from the H 1 -NMR spectra, where 
i— 1, 2, 3..., 200 (For raw data, see Supporting Information: table 
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Figure 1. 1H-NMR spectra. Purple lines represent the control group, red lines represent the mild injury group, green lines represent the severe 
injury group (3 days after injury), and blue lines represent the severe injury group (7 days after injury). Significant and distinguishable metabolomic 
changes occur after SCI. The right subplot shows the relative change in choline, lactate, lysine, and methylamine. 
doi:1 0.1 371 /journal, pone.0093736.g001 
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RMSE verified the validity of the PLS-regression model and was 
calculated according to the following formula: 



RMSE = \ - (observedi — predictedi) 2 

Vf=t 

RMSE = 4.76 indicated that the model had an acceptable 



Table 1. Summary of Eigen-metabolome: metabolites related 
to severe SCI. 





Categories 


Metabolite 


Fold-change 


P 


Neural Injury Metabolism 


4-Hydroxybutyric acid 


1.2 


ns 




Choline 


0.7 


<0.01 




L-Arginine 


1.2 


ns 




L-Histidine 


1.0 


ns 




L-Methionine 


2.1 


<0.01 




Tyramine 


2.3 


<0.01 




L-Serine 


1.8 


ns 


Energy Flow 


L-Lactic acid 


2.4 


<0.01 




Malonic acid 


1.2 


ns 


Skeleton Muscle 


Phosphocreatine 


0.8 


<0.01 




3-lndolebutyric acid 


0.7 


<0.01 




3-Methylhistamine 


1.3 


ns 




Putrescine 


2.1 


<0.01 


Urine Impairment 


Urea 


2.1 


<0.01 




Inosine 


1.5 


0.03 



Abbreviation: ns = no significance. 
doi:10.1371/journal.pone.0093736.t001 



accuracy (for raw data, see Supporting Information: table S2.) 
Following this, the relationship between the discrimination 
equations was used to evaluate the relationship between plasma 
metabolites and injury severity, where x represented NMR spectra 
ppm values and a t j represented loadings. 
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Finally, we arrived at an injury severity discrimination model 
based on the support vector machine (SVM). 



PC 2 + 1.075PCi- 1.269 >0 
PC 2 + 1.091PCi-1.114<0 
- 1.075PQ + 1.269 <PC 2 < 1.091PQ 



(mild injury) 
(severe injury) 
■1.114 (control) 



In summary, plasma NMR spectra could be clustered into one of 
three categories: control, mild injury, and severe injury. Only one 
mild injury sample was incorrectly clustered into the control 
category. All severely injured rats were separated from control and 
mildly injured rats (Figure 4). 

4. Metabolic network disturbance visualization after SCI: 
KEGG mapping 

After metabolites related to severe SCI were identified, 
metabolite roles and interactions were investigated on the global 
metabolic network. A three-step method was developed and 
included the following: 1) HMDB-based identification of enzymes 
that catalyze metabolites; 2) enzyme ID conversion from HMDB 
to KEGG enzyme IDs; 3) the use of Matlab KEGG API to map 
the enzymes into metabolic pathways (KEGG metabolic pathway 
ID: PATHWAY: mapOHOO). Following this, disturbed metabolic 
pathways and metabolic pathways related to neurobehavioral 
recovery after severe SCI were visualized (Figure 5). 

Seven metabolic pathways related to severe SCI and eight 
metabolic pathways related to neurobehavioral recovery were 
identified (Table 3). 
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Figure 2. Metabolic disturbance visualization after SCI: study design. (A) 1 H-NMR spectra (600 MHz) of plasma treated with a piecewise 
integral to form a 29x200 matrix. (B) PCA and PLS. Variable importance for the projection (VIP) values were used to screen metabolites. (C) The 
Human Metabolome Database (HMDB) was used to identify enzymes and genes that metabolize and regulate key metabolites. (D) Mapped 
metabolites and their enzymes at a global metabolic pathway in KEGG formed the visualization of metabolic disturbances after SCI. 
doi:1 0.1 371 /journal. pone.0093736.g002 



Discussion 

We successfully use plasma 1 H-NMR spectra to construct an 
Eigen-metabolome-based quantitative SCI severity evaluation 
model, which poses clinical application potential in future. Fifteen 
metabolites made up an Eigen-metabolome that distinguished 
severe SCI from mild injuries and healthy controls, whereas 
sixteen metabolites made up an Eigen-metabolome that distin- 
guished neurobehavioral recovery rats from other treatment 
groups. 

We established a PLS linear formulation group to provide a 
severe SCI evaluation model. In our model, X represented the 
ppm value from 1 H-NMR spectra and Y represented the severity 
of SCI. We found that we could distinguish control, severe, and 
mild injury rats in the PLS score plot by using SVM. To avoid the 
challenge by excessive clustering and small sample size, we used 
RMSE to evaluate the accuracy of our model. RMSE of our 
model (4.76) indicated it was appropriate for evaluating injury 
severity. 

Although few studies demonstrate that SCI disturbs plasma 
metabolites and metabolic patterns, SCI metabolic disturbances 
continue to challenge researchers. Consistent with our study, 
Fujieda Y et al. found that neurobehavioral recovery after SCI is 
related to a decrease in N-acetyl-aspartyl-glutamate and N-acetyl- 
aspartate in a liquid chromatography-mass spectrometry (LC— 
MS) metabolomic model [35]. However, this study focused their 
analysis on a few biomarker(s). To approach comprehensive 



understanding the complex network after SCI, to identify a few 
biomarkers is not enough. 

We introduced a new approach that integrates metabolomic 
analysis with metabolic network analysis. Significantly altered 
metabolites were extracted from the Eigen-metabolome and 
mapped onto a metabolic network. Enzymes that produce these 
metabolites were also identified and included in the network. 
Hundreds of plasma metabolites are altered after major trauma 
because enzyme activity increases or decreases to effectively 
decrease or increase metabolites levels [18]. We found that the 
following metabolic paths were disturbed: 

1) Retinal metabolism (Figure 5B) 

The retinal metabolism metabolic pathway contains two 
important metabolites, retinaldehyde (vitamin A) and retinoic 
acid (RA). RA plays a role in neural differentiation during 
embryonic differentiation. Studies have shown that RA possesses a 
neuroprotective effect during central nervous system injury. 
Plasma vitamin A, vitamin E, and beta-carotene are decreased 
after SCI. Furthermore, a positive correlation between plasma 
vitamin A and the severity of the chronic SCI has been identified 
[22]. A previous study found that RA at the injury site is increased 
after SCI [23]. 

2) Glycerophospholipid and arachidonic acid (AA) 
metabolism (Figure 5D-E) 

Phosphoglyceride and sphingomyelin could be used to distin- 
guish the severity of injury in 1 H-NMR spectra. Both phospho- 
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Table 2. Summary of Eigen-metabolome: metabolites related 
to neurobehavioral recovery. 





Categories 


Metabolites 


Fold-change 


P 


Inflammation 


Trimethylamine 


1 .6 


<0.01 




3-Methyhistamine 


1.3 


ns 


Glucose Metabolism 


Galactonic acid 


1.8 


<0.01 




Gluconic acid 


1.3 


ns 




D-Xylulose 


0.5 


<0.01 


Neurotransmitters 


Theophylline 


1.1 


ns 




Dopamine 


0.8 


<0.01 


Lipids 


Heptadecanoic acid 


1.3 


0.0318 




Pentadecanoic acid 


1.1 


0.047 




Glyceraldehyde 


2.1 


<0.01 


Amino Acids 


Dimethyl glycine 


0.9 


ns 




L-Serine 


1.1 


ns 




L-Glutamic acid 


1.2 


0.05 




L-Histidine 


1.0 


ns 




Oxalacetic acid 


0.9 


ns 




5-Aminolevulinic 


1.3 


<0.01 



Abbreviation: ns = no significance. 
doi:1 0.1 371/journal.pone.0093736.t002 



lipids are regulated by superoxide dismutase (SOD) and glutathi- 
one (GSH). SCI is currently considered as a myelin sheath 
damage. The myelin sheath is a greatly extended and modified cell 



membrane that is wrapped around the nerve axon in a spiral 
manner. Sphingomyelin and glycerol phospholipid are important 
components of the cell membrane, and when their metabolic paths 
are disturbed, the cell membrane is damaged. The simple 
sphingolipid metabolites ceramide and sphingosine- 1 -phosphate 
are major signaling cascade mediators involved in apoptosis, 
proliferation, stress responses, necrosis, inflammation, autophagy, 
senescence, and differentiation [19,20]. Consequently, cell mem- 
brane damage releases free fatty acid, particularly AA, from cell 
membranes. AA leads to secondary spinal cord neuron damage. 
The results of our study strongly support the hypothesis that free 
fatty acids contribute to tissue injury following spinal cord trauma 
[21]. 

3) The NAD-NADPH conversion process (Figure 5C, E, 
and I) 

The NAD-NADPH conversion process includes etheral 
metabolism, ethanol-ethylal conversion, and choline to betaine 
aldehyde metabolism. Many anabolic pathways such as lipid 
synthesis, cholesterol synthesis, cellular respiration via the electron- 
transport chain, and fatty acid chain elongation are related to 
NAD and NAD PH. Metabolomic fingerprints change after SCI, 
partly because of alterations in dehydrogenase activity. 

4) Tyrosine metabolism (Figure 5G) 

We found that the disturbances in tyrosine metabolism 
contributed to injury severity and neurobehavioral recovery. 
Tyrosine metabolism involves various neurotransmitter synthesis 
pathways. Tyrosine can be converted to L-DOPA, which is 
converted to the nerve injury biomarkers dopamine, norepineph- 
rine (noradrenaline), and epinephrine (adrenaline). The produc- 
tion and elimination of nitrotyrosine (NT) are markers of nerve 
injury during tyrosine metabolism. NT is a product of tyrosine 
nitration mediated by reactive nitrogen species such as the 



Variables/Loadings Plot 

PLS-DA Samples/Scores Plot 




Y Measured 1 



D E 

Figure 3. 1 H-NMR-based metabolomics spectrum: PCA and PLS. Pattern recognition of samples scored by PCA (A), PLS-DA (B and C), and PLS 

regression (D and E). (A) The PCA score plot of 1H-NMR spectra from plasma of controls (▼), severe injury (*), and mild injury (+). The ellipse was 
defined as the Mahalanobis distance from the distribution of samples. (B) Cross-validated scores plot of PLS-DA of the plasma 1H-NMR spectra for 
control (▼), mild injury (+), and severe injury (■). (C) VIP plot of PLS-DA. The threshold line indicates VIP value >1. (D) Regression plot of PLS-R based 
on BBB scoring of animals before scarification; the red line indicates the best fit line, the straight green line represents a 1:1 ratio, and the black line 
represents a 95% confidence interval. (E) Loading plot of PLS-R based on BBB scoring. 
doi:1 0.1 371 /journal. pone.0093736.g003 
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Figure 4. Severe and mild SCI score plot. The blue ellipse illustrates 95% CIs. The red and green lines represent SVM discriminant values. 
Multiclassification of severe injury (*), mild injury (+), and control samples (A). The SVM linear discriminant surface is composed of the red and green 
lines. LV = latent variables. 
doi:1 0.1 371 /journal. pone.0093736.g004 



peroxynitrite anion and nitrogen dioxide. NT is one of the most 
important markers for evaluating tissue and cell damage. Tyrosine 
nitration leads to enzyme inactivation [22,23] and affects tyrosine 
phosphorylation regulation, an important pathway in intracellular 



signal transduction [24,25]. Tyrosine nitration also plays a role in 
apoptosis [26] . Many diseases such as ischemia-reperfusion injury 
[27], chronic rejection after organ transplantation [28], cardio- 
vascular complications of diabetes [29], and Alzheimer's disease 




Figure 5. Metabolic networks disturbed by SCI: KEGG mapping. Nodes on the KEGG metabolic pathway for SCI severity (red) and BBB scoring 
(green) are shown as a series of disturbed metabolic modules in the network. (A) Estrogen and androgen transformation; (B) retinal metabolism; (C) 
enteral metabolism; (D) glycosphingolipid biosynthesis; (E) glycerophospholipid metabolism; (F) pentose and glucuronate interconvertion, 
glycoaminoglycan biosynthesis; (G) tyrosine and caffeine metabolism (H) choline to betaine aldehyde, serotonin metabolism (I); ethanol to ethanol, 
chloroalkene and chloroalkene degradation; and (J) cadaverine and putrescine metabolism. 
doi:1 0.1 371 /journal. pone.0093736.g005 
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Table 3. Metabolic pathways related to severe SCI and neurobehavioral recovery. 



Metabolic Pathways Related to Severe SCI and Neurobehavioral Recovery 





Severe SCI 


Neurobehavioral Recovery 


Shared Metabolic Pathways 


1. Glycerophospholipid metabolism 






2. Etheral metabolism 






3. Tyrosine metabolism 




Distinct Metabolic Pathways 


1. Retinal metabolism 

2. Chloroalkene and chloroalkene degradation 

3. Choline to betaine aldehyde metabolism 


1. Galactose and 4,6-phospho-D-gluconate metabolism 

2. Cadaverine and putrescine metabolism 

3. Caffeine metabolism 

4. Histamine to histidine metabolism 

5. Taurine and hypotaurine metabolism (Tryptophan metabolism) 



Note: Three metabolic pathways are related to both severe SCI and neurobehavioral recovery. 
doi:1 0.1 371/journal.pone.0093736.t003 



[30] are accompanied by high NT levels, further highlighting the 
importance of protein tyrosine nitration in tissue damage. 

5) Cadaverine and putrescine metabolism (Figure 5J) 

The role of decarboxylase ornithine decarboxylase (ODC) as a 
catalytic product in SCI has gained interest. ODC and putrescine 
levels significantly increase after SCI [3,31,32]. We hypothesize 
that coordinating these processes and other sequential processes 
such as spermine expression, enhancement of nitric oxide, and 
purinergic pathways will be effective to influence innate immune 
processes after SCI [33]. 

6) Disturbances in ether and alcohol metabolism 
(Figure 5C, H, and J) 

We found that several metabolic pathways such as chloroalkene 
and chloroalkene degradation and aldehyde and ethanol trans- 
formation (not well studied in SCI) appeared in the disturbance 
pathways. These metabolites belong to etheral metabolism and are 
primarily generated by the gut microbiome. This disturbance 
certainly takes place in clinical practice; however, it is typically 
ignored in traditional trauma biomedical research. 

In summary, SCI can be understood as disturbances in cellular 
signaling regulatory networks in response to trauma. Furthermore, 
these disturbances commonly present a pattern of general cellular 
transduction, as illustrated in Figure 3. It is possible to monitor and 
assess SCI cellular signaling and modulation changes by profiling 
the global metabolome and visualizing metabolic and gene paths 
into distinguishable patterns. 

Our study demonstrated that 1H-NMR spectra can be used to 
quantitatively measure and visualize the metabolic network 
disturbances after SCI. We identified major metabolites and 
enzymes through unbiased screening of compounds identified by 
metabolomic analysis. Metabolites were mapped to a metabolic 
network, and their functions were interpreted by KEGG metabolic 
network visualization. This method has the advantage of 
integrating metabolic research with the trauma mechanisms to 
aid our understanding of the pathophysiological SCI process. 

Our model incorporated metabolomics, biological informatics, 
and pattern recognition to develop a new evaluation and outcome 
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